Particle motion around magnetized black holes: Preston-Poisson space-time. 
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We analyze motion of massless and massive particles around black holes immersed in an asymp- 
totically uniform magnetic field and surrounded by some mechanical structure, which provides the 
magnetic field. The space-time is described by Preston-Poisson metric, which is the generalization 
of the well-known Ernst metric with a new parameter, tidal force, characterizing the surrounding 
structure. The Hamilton-Jacobi equations allow separation of variables in the equatorial plane. 
The presence of tidal force from surroundings considerably changes parameters of the test particle 
motion: it increases the radius of circular orbits of particles, increases the binding energy of mas- 
sive particles going from a given circular orbits to the innermost stable orbit near black hole. In 
addition, it increases the distance of minimal approach, time delay and bending angle for a ray of 
light propagating near black hole. 

PACS numbers: 04.30.Nk,04.50.+h 



I. 1. INTRODUCTION. 

Black holes in the centres of galaxies are immersed in 
a strong magnetic field due to charged matter surround- 
ing them. The strong magnetic field in the centre of 
galaxies is stipulated by toroidal currents around galactic 
black holes [l[ . Therefore an exact solution of Einstein- 
Maxwell equations describing a black hole immersed in an 
asymptotically uniform magnetic field, known as Ernst 
solution [U, was of considerable interest 0. The light 
and particle motion around Ernst-Schwarzschild black 
hole was analyzed in a few papers ||, In particular, 
in 0] it was shown that Hamilton- Jacoby equation allow 
separation of variables in the equatorial plane, where the 
motion of neutral and charged particles were analyzed. 
In Q the motion of neutral particles were considered for 
a more general situation of electromagnetized Kerr back- 
ground. There it was shown that the release of binding 
energy is considerably increased because of presence of 
electromagnetic field, and the binding energy for circular 
orbits was calculated. Yet, in a more realistic situation, 
the strong magnetic field in the central region near black 
hole is created by some surrounding matter, such as ac- 
cretion disk or an active galactic nuclei. This surrounding 
structure exerts strong gravitational tidal force on parti- 
cles moving near black holes, so that magnetic influence 
of the structure might be even much smaller than its 
gravitational influence. Therefore a more physical situa- 
tion should include into consideration, the corrections to 
the black hole metric due-to that structure. Fortunately, 
recently Preston and Poisson Q have found such a cor- 
rected metric. This is the solution to the perturbative 
Einstein-Maxwell equations depending on three parame- 
ters: the black hole mass M, magnetic field B, and a new 
parameter K, which characterize the above surrounding 
structure. The solution is very accurate for r 2 B 2 <C, 



M/a <C 1, and r 2 K <C 1, where r is the distance form 
the black hole, a is the length scale of the mechanical 
structure. Indeed, comparison with the exact Ernst so- 
lution shows that next order corrections are of order i? 4 , 
and are very small. 

In thepresent paper we generalize analysis of works 
[f| and Q and study motion of test particles near black 
holes immersed in asymptotically uniform magnetic field 
and some gravitating surrounding structure, which pro- 
vides the magnetic field. The paper is organized as fol- 
lows: In Sec. II we reduce the Preston-Poisson metric 
to the Ernst-like form, by going over to a new coordi- 
nates. Then in Sec. Ill we consider the Hamilton-Jacobi 
equation in the equatorial plane, and use it for analysis 
of massless particles. There the lens effects for Preston- 
Poisson metric is considered. Motion of massive particles 
is described in Sec. IV, where the binding energy for par- 
ticles on circular orbits are calculated. 



II. PRESTON-POISSON METRIC 

Following let us consider a model consisting of a 
non-rotating black hole immersed in a uniform magnetic 
field, and a large mechanical structure, such as a giant 
solenoid, producing the magnetic field of strength B. The 
structure has a mass M' and its linear extension is ~ a. 
In order to have magnetic field which is uniform for r S> 
M, one chooses 



A a ±B<f> a , 



(1) 



where <\> a is the rotational Killing vector of the unper- 
turbed Schwarzschild metric. 

The metric which describes the space-time of the above 
model, written in light-cone gauge, is [3| 
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ds 2 



-g vv dv 2 + 2dvdr + g v gdvd9 + gggdO 2 + g^dfi 2 , (2) 
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where 



g vv = -(1 - (2M/r)) - ^S 2 r(3r - 8M) 



(l/3)r 2 (3M -r)(M + r)B 2 , (12) 
= r 2 - (l/3)r 2 (-6M 2 + AMr + 5r 2 )K- 



B\3r z -lAMr + lMr) + K{r-2My ) (3cos 2 6»-l), 

(3) 



B 2 (r - Ml) -2K)r 2 sin 9 cos 9 (4) 



(l/3)f 2 (3M 2 - 2Mr + r 2 )B 2 



(13) 



This form of the metric does not have non-diagonal com- 
ponents in the equatorial plane and is much simpler 
for consideration of motion of test particles. To be ex- 
act, non-diogonal components are of order (B 4 , K 2 ) and 
higher, and therefore can be safely neglected. 



r 2 + ( -^B 2 r 2 + B 2 M 2 + K{r 2 - 2M 2 ) ) r 2 sin 2 ■ 



(5) 



94,4, = r 2 sin 2 9+ 



B 2 r 2 [y-\r 2 - M 2 j - K(r 2 - 2M 2 ))r 2 j sin 4 9 (6) 

The above metric is accurate through order (B 2 , K), 
whenever r 2 B 2 <C and r 2 K <C 1. The new parameter K 
characterizing the mechanical structure, containing the 
black hole, can be interpreted as a tidal gravity or Weyl- 
curvature when r S> M. Thus we have a three-parameter 
solution. 

From now on, we shall consider motion in the equato- 
rial plane 9 = tt/2. Therefore, we shall put 9 = ir/2 in 
formulas (l)-(5). Then, let us make the following coor- 
dinate transformations: 



v = t + 7+ 2M In I (r/2M) - 1|, 



(7) 



r = r(l + (l/6)B 2 r 2 + (1/3) (# - (l/2)£ 2 )(r - 2M)r 2 + 



0(B 4 ,K 2 ). 



(8) 



These transformations are generalization of transforma- 
tions (3.60-3.61) of [4| in the equatorial plane and they 
cast the metric (1) into the diagonal form 

ds = gttdt + grrdr + geed9 + g^dcj) , (9) 

where the metric components (neglecting orders 
0(B 4 ,K 2 ) and higher) are 

, , , (AM - 3r)r 2 2Mr 2 9 

opp = (1 - 2M r)~ 1 + - — K B 2 , 

y y ' ' 6M-3r 6M-3r 

(10) 



g tt = (1 - 2M/r) + (1/3) (-8A/ 2 + VQMr - 3r 2 )K+ 



(2/3)M(2M- r)B 2 , 



r 2 + (l/3)r 2 (-6M 2 - 4Mr + 5r 2 )K+ 



(11) 



III. MOTION OF MASSLESS PARTICLES 

From now and on we shall write r instead of r. The 
four-momentum is 



dx^ 
ds 



(14) 



where s is an invariant afhne parameter. The Hamilto- 
nian has the form 



H = \g^ P ^. 



(15) 



The action can be represented in the form: 

S= -fis- Et + L4> + S r (r) +S 6 {9), (16) 

where E and L are the particle's energy and angular mo- 
mentum respectively. 

Then, the Hamilton- Jacobi equations for null geodesies 
read 



1 av 9S dS 
— cr v 

2 dx» dx v 



dS_ 
ds 



(17) 



It is evident that the equations of motions allow sepa- 
ration of variables in the equatorial plane 9 = tt/2. The 
first integrals of motion are 



dr 
ds 



Eh 

gu 



Mrr~-=± X \-»—{E 2 ~Ul ss X 



P<t> = = L , 



dt 

Pt = H9U-T- = ~E, 
as 



u; ff = -gtt^ 1 + — 

M 9H 



(18) 

(19) 
(20) 
(21) 



The trajectory and propagation equations take the 
form 



gtt L gtt (PtY m2 
g rr \ g^ \pt ) u p 2 



(22) 
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fdrV _ g<M, ( 944 ( pt\ 2 „ ™ 2 \ (9 o^ 

For massless particles, from the above equations (22), 
(23), one can see that propagation and trajectory equa- 
tions contain only ratio b = L/E, which is called the 
impact parameter. The qualitative description of the mo- 
tion of massless particles can be made by considering the 



effective potential of the motion (21, where /i — 0. The 
equation for radii of circular orbits can be found from the 
condition dU e ff/dr = 0: 

gttg4>4>,r = g<j,<j,gu,r ■ (24) 

This gives the algebraic equation for r: 



6(B 2 - 2K)M 6 - 2r + 10(B 2 ~ 2K)M 2 r + (4/3) {B 2 + K)r + M(6 - 65V + AKr 



0. 



(25) 



When B 2 -C M, and K -C M, we have, that above some critical region of values of B and K, there are two null 
circular orbits with radii 



n = (3M + 12M 3 ,ft: + (3M 3 + 120M 5 K)B 2 



(26) 



r-2 



2K(y/G - 6y/KM)(l + 2KM 2 ) - B 2 (V6 + y/KM(- 15 + 13y/6KM + 6KM 2 )) 

A{Kf/ 2 ' 
I 



(27) 



When B = K = 0, r*i takes its Schwarzschild value 
3A/. Unfortunately, we cannot find accurate values for 
the critical region (B cr , K cr ), because the values we get 
is quite large, and, for instance for K = B 2 /2 is about 
0.189M, what is on the boundary of applicability of the 
approximate metric under consideration. Physical situ- 
ation corresponds to some tidal force K, which is larger 
than its pure Ernst value B 2 /2. Therefore, further we 
shall consider the new parameter h, which is given by 
the relation: 



K = 



B 2 



h. 



Now, let us consider the effect of tidal gravitational at- 
traction of the surrounding structure upon such lens ef- 
fects as light bending angle and time delay. For this, let 
us follow the approach of 0- 

If we know the distance of minimal approach rw n with 
great accuracy, we can perform integrations for finding 
bending angle: 



—dr 

dr 



d0 

—dr — it. 

dr 



(28) 



Here r a is radial coordinate of an observer and r s is 
radial coordinate of the source. 

In a similar fashion one can find the time delay, which is 
the difference between the light travel time for the actual 
ray, and the travel time for the ray the light would have 



Table I: Bending angle a and propagation time r for Ernst- 
Schwarzschild space-times (in geometrical units, M — 1) for 
6 = 6. "Observer" and "source" are supposed to be situated 
not far from the black hole in order to estimate influence of 
magnetic field in the central region of the black hole: r = 
r s = 20, b = 6. 



B 


h 


a + 7r 


ts—to-b d s -o/ cosS 








4.252334 


56.84554 





5- 10" 4 


4.406231 


57.25224 





10~ 3 


4.577176 


57.72044 


5- 10~ 4 





4.252388 


56.84567 


5- 10" 4 


5- 10" 4 


4.406293 


57.25240 


5- 10" 4 


10~ 3 


4.577247 


57.72063 



taken in the Minkowskian space-time: 



ts to — 



dt 
dr 



dr 



— dr - ds ~° 
dr cos B ' 



(29) 



Here the term c ^~^ represents the propagation time for 
a ray of light, if the black hole is absent. The distance of 
minimal approach is the corresponding root of the equa- 
tion dr/dt = 0. For pure Schwarzschild black hole it 
would be the largest root, yet in our case it the second 
largest root of the equation: 

r 3 (-3 + B 2 (3M 2 - 2 Air + r 2 ) + K(-6M 2 + AMr + r 2 )) 

- (2M - r)(3 + 2B 2 Mr + Kr(—AM + 3r)b 2 = 0. (30) 

Looking numerically for the solution of the equation (30) 
one can see that the tidal force K pull out the radius 
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Figure 1: Effective potential for neutral particles: M = 1, 
B = 0.001, fi = 0.1, L = 0. K = 0.001 (blue), K = 0.01 
(green), K = 0.02 (red). 
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Figure 2: Effective potential for neutral particles: M = 1, 
£ = 0.001, a* = 0.1, L = 5. if = 0.001 (blue), K = 0.01 
(green), K = 0.02 (red). 



of minimal approach further from the black hole. We 
also can see from the table I, that the presence of the 
mechanical structure leads to increasing of the banding 
angle and time delay near the black hole. 



IV. MOTION OF MASSIVE PARTICLES 

The effective potential for massive neutral particles 
(21) is shown in Figures 1 and 2 for zero and non-zero 
angular momentum L of the particle. 

From the above figures one can see that the effective 



potential can have form of the barrier or of a monotoni- 
cally increasing function. 

For circular orbits the equation V e //, r = gives 



L\gt 



0. (31) 



If one solves (27) for L, and uses it in the equation 
dr/dt — 0, or, equivalently, in 



(32) 

one obtains rather cumbersome system of equation for 
determination of the parameters of orbits of massive par- 
ticles. 

The values L//J, and Ej \x as a function of radius of 
circular orbits are presented on Fig. 3 and 4 . They 
are found there from accurate equations (31) and (32). 
From Fig. 3 and 4, one can see that particle angular 
momentum and energy per units mass is monotonically 
growing as functions of the radius of circular orbit r c , 
starting from some minimal value. This minimum value 
of the test particle angular momentum corresponds to the 
orbit with the innermost stable circular radius ri C . The 
large K, the more radius of the innermost stable circular 
orbit ri c is pulled toward the black hole. 

The binding energy is defined as the amount of energy 
that is released by the test particle going from a stable 
circular orbit r c , to the innermost stable orbit of radius 
r ic . i.e. 



Binding energy 



{E/n) rc - {E/u) r 
(£/Mk 



(33) 



A test particle in an unstable circular orbits will fall 
into black hole and the infall time is small compared to 
the radiative time, so that the particle energy will be 
brought to the black hole almost completely. 

From the Figure (5) one can see the two features: first, 
that the binding energy is greater for larger radius of 
circular orbit, and this dependence on radius is strictly 
monotonic. Second, is that the larger tidal force K, the 
larger the binding energy for a given radius of the circular 
orbit of the particle r c . The last feature means that in 
presence of surrounding attracting structure a test par- 
ticle, when going from its stable orbit to the innermost 
stable one, releases more energy, than it would release 
without the above structure. 

Finally, let us give the explicit form of the equations 
(31), (32) through the orders B 2 , K. 



E 2 1 f r — 2M\ 2 /3 

" - ' > ' -(3M — r) + 2K(9M 2 - lQMr + 3r 2 ) — M(9M — 4r)B 2 I (34) 



^2 ~ 3 {r-3M 



[i 2 3 \r-3MJ \r 

2 l , , ^ 

(3r(r - 3M) - 18M 4 + 6M 3 r + l9M 2 r 2 - UMr 3 + 3r 4 + M(9M 3 - 3M 2 r - 2Mr 2 + r 3 )B 2 ). (35) 
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Figure 3: as a function of radius of the circular orbit r c 

M = 1, B = 0.0001, K = {B 2 /2) (bottom), K = (B 2 /2) + 
0.001, K = (B 2 /2) +0.003 (top). 




Figure 4: L/^t as a function of radius of the circular orbit r c 
M = 1, B = 0.0001, K = {B 2 /2) (bottom), if = (B 2 /2) + 
0.001, A" = (B 2 /2) +0.003 (top). 



Binding Energy 




Figure 5: Binding energy as a function of radius of the circular 
orbit r c M = 1, B = 0.0001, K = (B 2 /2) (bottom), if = 
(B 2 /2) + 0.001, K = (£ 2 /2) + 0.002, K = (B 2 /2) + 0.003 
(top). 
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Figure 6: Effective potential for charged particles: e = 0.6 
(top), 0, -0.6, -3.8 (bottom), h = 0, fiM = 1, L = 30 



These expressions are much simpler than exact equa- 
tions (31), (32). 

For massive charged particles in vicinity of uncharged 
black holes, one should change the angular momentum 
of the particle L, to the generalized momentum: 

L -> L + g^^—. (36) 

The effective potentials for the case of charged particles 
are shown on Fig. (6)-(8). The situation is dependent on 
the sign of the charge, because the Lorentz force acting on 
the charged particles has opposite directions for positive 
and negative charges. 

There are two reasons why we do not analyze the case 
of charged particles in detail: First, the magnetic field 
used for derivation of the considered metric is given only 
trough the first order in B. Second, the effect of strong 
magnetic field for charged particles is stipulated by the 
factor eB/n (when M = 1), and is very large even for 
small B -C M, because of the large ratio e/fj,. Therefore, 
it is generally accepted, to neglect "geometric" influence 



on propagation of charged particles, and to consider the 
more realistic decaying magnetic fields on the black hole 
background @|. 

V. CONCLUSION 

We have considered the motion of massless and mas- 
sive test particles near black holes immersed in asymp- 
totically uniform magnetic field and some surrounding 
structure which provides this field. The tidal force from 
the surrounding structure has considerable influence on 
the parameters of the test particle motion. Let us enu- 
merate them: a) it pulls radius of the circular orbits off 
the black hole, b) increases the radius of minimal ap- 
proach for light, c) increases the time delay and bending 
angle for light, d) increases the energy and momentum 
(per unit mass) for a circular orbit of a given radius, e) 
increases the binding energy of massive particles, which 
releases when a particle goes from a given stable circular 
orbit to the innermost stable circular orbit, f) the radius 
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of the innermost stable circular orbit is pulled closer to 
the black hole. 

The used Preston-Poisson metric gives an excellent op- 
portunity to investigate the motion of test particles in the 
vicinity of a supermassive "dirty" black hole, surrounded 
by some distribution of matter and uniform magnetic 
field, and to approach, thereby, a more realistic situa- 
tion than that given by the Ernst solution. 



Figure 7: Effective potential for charged particles: e = 0.6 

(top), 0, -0.6, -3.8 (bottom), h = 0.005, fiM = 1, L = 30 U eff 




Figure 8: Effective potential for charged particles: eB = 0.006 
(top), 0, -0.006, -0.038 (bottom), h = 0, \iM = 1, L = 5 
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